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The Mesoscopic Mechanics (MeM), as introduced in [5], is relevant to the 
electron gas confined to two spatial dimensions. It predicts a special way 
of collective response of correlated electrons to the external magnetic field. 
The dynamic variable of this theory is a finite-dimensional operator, which 
is required to satisfy the mesoscopic Schrodinger equation, cf. (2) below. 

In this article, we describe general solutions of the mesoscopic Schrodinger 
equation. Our approach is specific to the problem at hand. It relies on the 
unique structure of the equation and makes no reference to any other tech- 
niques, with the exception of the geometry of unitary groups. In conclusion, 
a surprising fact comes to light. Namely, the mesoscopic dynamics "filters" 
through the (microscopic) Schrodinger dynamics as the latter turns out to 
be a clearly separable part, in fact an autonomous factor, of the evolution. 
This is a desirable result also from the physical standpoint. 

1 A brief description of contents and results 

The mesoscopic Schrodinger equation describes evolution of an operator (de- 
noted K) via a nonlinear equation. In order to motivate the reader let me 
point out that the interpretation of this operator is somewhat similar to 
that of the wavefunctions (of the regular Schrodinger equation). Informally 
speaking, the pair of operators K and K* may be interpreted as essentially 
being square roots of a density matrix, (cf. [5]), although the issue is del- 
icate due to non-commutativity. This should be viewed as comparable to 
the fact that modulus-square of a wavefunction represents a probability dis- 
tribution. I have proposed the mesoscopic Schrodinger equation as a model 
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for correlated evolution of an n-tuple of electrons, which is relevant to the 
galvanomagnetic properties of the so-called correlated materials. 

In this article, I present analysis of the evolution of the system when the 
single particle Hamiltonian, which is an ingredient in the equations, does not 
receive any feedback from the dynamic variable K. The last section of this 
article, Section 5, provides a brief summary of the physical interpretation 
of the mesoscopic equation. In particular, it should explain why there is 
incentive also to consider the case when such a feedback would exist. This 
problem is not addressed in the present article. However, as regards the 
case limited to the i^-independent Hamiltonian H, the problem is essentially 
resolved. 

Section 2 is meant to introduce the system, and to review some basic 
properties, display simple special solutions. Also, it is pointed out that the 
mesoscopic equation has the structure of a Hamiltonian system. However, I 
emphasize, no further use is made of the so-called canonical formalism. Next, 
in Section 3, the equation is solved in the case when the domain and image of 
the operator K are fixed finite-dimensional spaces. In fact, it is shown that 
the evolution of this nonlinear system with time-dependent Hamiltonian can 
be represented in a certain way via a pair of curves on the unitary group. 
Finally, we consider the case when the domain and the image of operator K 
are a priori allowed to evolve in an ambient Hilbert space. In this case, the 
single-particle Hamiltonian is densely defined on the Hilbert space. Solutions 
of the mesoscopic equation in such a broad setting are described in Section 4. 
In particular, a uniqueness property is shown. Its proof takes advantage of an 
exceptional structure of the mesoscopic equation and it could not, it seems, 
be deduced from any general principles. Also, it is shown that the evolution 
is "driven" by an n-tuple of Schrodinger particles. This is important from 
the physical stand-point as obtaining any other type of carriers might be 
problematic from the standpoint of physical interpretation. In addition, we 
note that evolution entails a phase factor which explicitly depends on the 
history of the magnetic energy density B 2 (t). Thus, in the fixed-domain and 
the moving-domain cases alike, solutions are represented by means of simpler 
factors corresponding to certain linear problems. We emphasize that this 
remains valid even when the constituents of the equation are time dependent. 
This is what is meant by integrability in the mesoscopic dynamics. 

Let me emphasize that while considerations in Section 3 are related to 
the ordinary differential equations, those of Section 4 deal with partial differ- 
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ential equations. However, I believe, the context of operator equations with 
evolving domains may be quite new. It enables one to capture essentially 
new phenomena that cannot be discussed on grounds of the PDE setting 
alone. 

2 The mesoscopic Schrodinger equation 

Let F be a finite-dimensional complex vector space equipped with a Hermi- 
tian scalar product. Let t denote the time variable, and let 

H(t) :F^F, 

be a predetermined family of positive definite self-adjoint operators, which we 
will refer to as the Hamiltonian. In addition, let B — B(t) be a predetermined 
function of time, which we will refer to as the magnetic induction. We 
will require throughout this article that both B and H depend on the time 
variable smoothly. This is a technical assumption, which will ensure local 
existence and uniqueness of solutions of certain dynamical systems that we 
will encounter along the way. Introduce the dynamic operator variable 

K(t) :F^G. (1) 

We assume that K has a null kernel, kerfT = {0}, while the target space 
G = lm(K(t)) is an arbitrary but fixed complex linear space also equipped 
with a Hermitian scalar product. (In fact, we will consider a more general 
situation in Section 4.) Throughout this article our attention is focused on 
the mesoscopic Schrodinger equation 

ihk = -KH - B 2 (K*)-\ (2) 

where the * denotes Hermitian conjugation. Note that the nonlinearity is of 
a homogeneous type but develops a singularity as det K — > 0, which a priori 
may be an intimidating factor as one attempts to solve the equation. 

Let us recall here that the manifold of invertible linear transformations, 
say, from F to G, is equipped with a natural Hermitian metric given by 

(L\N) = trace (LN*) . 
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Here L and N denote two arbitrary tangent vectors which, let it be empha- 
sized, represent arbitrary linear transformations from F to G. Furthermore, 
the Hermitian structure induces a compatible Riemannian structure 

(L,N) = 3? {trace (LA^*)}, 

as well as a symplectic form 

u(L,N) = 3 {trace (LN*)} . 

With this understood, let us point out that the evolution equation (2) is tied 
to the following total Hamiltonian 

E(K) = trace (KHK*) + B 2 log det (KK*) . (3) 

Indeed, a calculation shows that the differential of S is given by 

dZ K [L} = i\ £=Q Z(K + eL) = dt {trace ((KH + B 2 (K*)~ l ) L*)} 

= (KH + B 2 (K*)-\L), { ] 

Furthermore, since 

n {trace (AB*)} = 3 {trace (iAB*)} , 

equation (4) can be re-interpreted in the form 

dE K [L] =oo [i (KH + B^K*)- 1 ) , L) . 

This means precisely that (2) is the Hamiltonian flow (cf. [1]) induced by 
the total Hamiltonian S and the symplectic structure uj. While it is good to 
bring this theme to the reader's attention, it will not be explicitly essential 
to the discussion in this article. 

Formula (4) indicates that the critical points of the Hamiltonian S, sub- 
ject to the constraint trace {KK*) = const, satisfy the Euler-Lagrange equa- 
tion 

KH + B 2 (K*)~ l = vK. (5) 

This equation implies 

K*K{v -H)= 5 2 Id. 
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In addition, since K*K > 0, the equation can be satisfied only if the real 
scalar v dominates all the eigenvalues of H . In conclusion, all solutions of 
(5) are of the form 

K > = u J7^Wn' < 6) 

where U : F — > G is an arbitrary unitary operator, and v is arbitrary as long 
as it dominates H. The critical points are interesting in their own right, cf. 
[5]. In addition, they play a special role in the time-dependent problem (2). 
Indeed, assume for a while that the Hamiltonian H is time-independent and 
diagonalized by vectors \tp n ) G F, so that 

H\^ n )=E n \^ n ) (7) 

for a collection of positive eigenvalues E n . Apparently, the simplest solutions 
of equation (2) are of the form 

K = J2an(tMn)(^n\, (8) 

where a n = r n e l(fn . Substituting this into (2) one readily obtains 



I / B 2 \ 

r„ i0 and ip n = - I E n + — )t + (p nfl . (9) 

II V r n,0/ 



It is interesting to note that when all correlated, i.e. oscillate with 

the common frequency, say, v = E n + B 2 jr\ for all n, then 



±B 

r n = 



[y - En) 1 / 2 ' 



Therefore, these special solutions conform with (6), and so they represent 
critical points of the Hamiltonian S. 



3 General fixed- domain solution of the meso- 
scopic equation 

In this section, it will be shown that equation (2) can be reduced to a sys- 
tem of simpler equations, even when H and B, i.e. the constituents of the 
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equation, are time dependent. As we set out to solve the equation, the first 
useful artifice is to use polar representation of the operator. Namely, let 



K = RU, 



(10) 



where R = R* is positive definite, and U~ x = U*, i.e. U is unitary. It ought 
to be emphasized that here the matrix U stands on the right, which is in 
contrast to the situation in (6). For a given K, its polar representation is 
determined by setting 

R = VKK* : G -> G, and U = R~ l K : F -> G. 

A direct calculation shows that U selected in this way is unitary. It is well 
known that with the requirement of positive definiteness of R the polar de- 
composition is unique. Next, observe that when K satisfies equation (2), 
then 



%KK* + kKK* 

—i (-KH — B 2 (K*y l ) K* 



iKHK* + iB 2 
0. 



iKHK* 



+ tK( 
iB 2 



-HK* - B 2 K- 1 ) 



It ought to be emphasized again that the calculation remains valid whether 
or not the Hamiltonian and the magnetic field depend on the time variable. 
Let us now set K(0) = K . We have 

R 2 = K K*. 

There is only one positive definite, self-adjoint R satisfying this condition. 
Since K Kq is positive definite and self-adjoint, it can be diagonalized in a 
certain basis so that 



tf tf * = diag XlXl 



\ 2 



and, in the same basis, 



i? = diag[|Ai|,|A 2 |,...,|Ajv|]. 



R does not depend on time. In summary, we obtain 
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Corollary 3.1 Evolution prescribed by equation (2) is constrained to the 
submanifold 

M R = [K : KK* = R 2 }. 

All coordinate functions of the matrix KK* are integrals of motion. As is 
easily seen, M R is diffeomorphic with the unitary group and has half the 
dimension of the phase space. 

A similar calculation as above shows that 

^-E(K(t)) = trace (KHK*) + 2BB log det R 2 . (11) 

In particular, as the system evolves, change in the entropy part of the total 
Hamiltonian only depends on B(t). If H does not depend on time, then 
trace [KHK*) is an additional integral of motion. 

We now continue to discuss solutions of (2). First, denote 

17(0) = R~ l K = U . 
Next, substitute K in its polar representation into equation (2) to obtain 

ihRU = -RUH - B 2 R- l U 
Multiplying the equation by R~ l one further obtains 

ihU = -UH - H B U, (12) 

where 

H B = B 2 R- 2 = B 2 (K Q K%)- X . (13) 

In this way, evolution of the unitary part is determined by the predetermined 
constituents H and B as well as the initial condition K . In fact, it may be 
more practical for some purposes to represent equation (12) in the form 

U*U = Uh + U*H b U). (14) 

The left-hand side represents a vector tangent to the trajectory, shifted to the 
group unit. The right hand side, driving the evolution, represents an element 
in the Lie algebra of skew-Hermitian operators. Indeed, not only H but also 
H B , and hence also U*HbU are Hermitian operators. Suppose at first that 



7 



H and B are all frozen in time. Since the unitary group is compact and the 
group multiplication is smooth, the right-hand side of (12) defines a Lipschitz 
continuous vector field on the unitary group. In particular, solutions of (12) 
are uniquely determined (via a choice of the initial condition) and exist for 
all time. In fact, in this case the solution may be written in the form of a 
power series 

U(t) = U + i {{U Q H + H B U ) - ^i(U H 2 + 2H B U H + H%U ) 

(15) 

-i^{U H 3 + 3H B U H 2 + 3H B U H + H B U ) + ..., 

Since all operators are finite dimensional, the series converges absolutely. A 
straightforward calculation shows that U(t) satisfies (12). 

We proceed to resolving the case when H and B are allowed to vary in 
time smoothly. First, represent U as a product of two unitary matrices, i.e. 

U(t) = V(t)W(t). 

Equation (12) yields 

ihVW + ihVW = - VWH - H B VW. (16) 

Secondly, multiply the equation by V* on the left and by W* on the right. 
This leads to 

ihV*V + ihWW* = -WHW* - V*H B V. (17) 

Now, the two factors have been separated. Indeed, ask that W and V satisfy 
the following two separate equations 

ihW = -WE, (18) 

and 

ihV = -H B V. (19) 

In the case of time- varying H and H B only local existence of solutions of (12), 
(18), and (19) is guarantied, but the uniqueness property is still retained. It 
follows that if 17(0) = V(0)W(0), then U(t) = V(t)W(t) for all t. Note that 
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Hb depends on time only via B, and due to Hermicity, it can be written in 
a certain basis as 



H B (t) = B 2 (t)dmg[\i 2 ,X^,...,X N 2 
Thus, the solution of (19) can be represented in the same basis in the form 
V(t) =exp^J^B 2 (t')(K K*y 1 dt') 

= diag [exp (^Ar 2 Jo* B 2 (t')dt>) , . . . , exp ^X~ N 2 J* B\t')dt>) . 

Here, we have selected the initial condition V(0) = Id. This needs to be 
compensated by the appropriate choice of the second initial condition, namely 
W(0) = Uq. As it turns out, we have essentially reduced equation (2) to a 
pair of simpler, well-understood equations. Let us summarize the results. 

Theorem 3.1 Consider the mesoscopic Schrddinger equation (2) with smooth 
constituents H = H{t) and B = B{t). The solution K = K{t) satisfying the 
initial condition 

K(0) = K = RU 

is a uniquely defined smooth operator-valued function of time. Furthermore, 
the solution admits representation in the from 

K{t) = ^jK K* exp ^-j^B 2 {t'){K Kl)- l dt^ W(t), (20) 

where W satisfies 

iKW = -WH(t), W(0) = U . 

When H and B are time-independent the solution exists for all time, while 
in general it is only guarantied to exist locally. 

Of course, if both Hb and H are time-independent, then V(t) and W(t) 
represent two geodesies of the bi-invariant metric on the unitary group, e.g. 
cf. [4]. They can also be represented as power series. One can perform 
multiplication of the two series and grouping of the terms to see that the 
product is equivalent to the series in equation (15). It is worthwhile to 
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mention that when H depends on time, W(t) can still be represented in 
terms of the time-ordered exponential, cf. [3] p. 219. 

It is worthwhile to substitute the solution of (2) in the form specified in 
Theorem 3.1 into formula (11). A calculation involving the property that 
trace (AB) = trace (BA) shows that the following holds. 

Corollary 3.2 In the notation of Theorem 3.1, we have 

^rZ(K) = trace (W*R 2 WH) + ^(B 2 ) logdet R 2 . (21) 

(JjL (JjL 

Since R is fixed in time, the magnetic (entropy) part of the energy only 
depends on magnetic induction during the evolution. Recall that evolution 
of W only depends on H , and so the electronic part of the energy is only 
affected by the electronic constituent. ( Of course, H could depend on B via, 
say, Landau quantization.) 

I would also like to highlight the fact that we have made many arbitrary 
choices when solving equation (2). Naturally, we have made those choices 
so as to simplify the discussion. In spite of that, uniqueness of solutions 
guaranties that the result is general. One of the very conspicuous arbitrary 
choices was declaring time-independent operator R. We need not impose the 
condition of positive definiteness of R. If that condition is dropped and when 
K Kq has degenerate eigenvalues, one can select a time- varying R satisfying 
the constraint R 2 = KqKq. Subsequently, one would redefine the auxiliary 
Hamiltonian Hb by setting Hb = B 2 R~ 2 + iUR^R. Naturally, this would 
also redefine V = V(t) and in the end yield the same product RV as the 
calculation based on the time-independent R. 

Let us look back at the findings in this section. Recall that equation 
(2) has a strong yet homogeneous nonlinearity. In fact, one might argue it 
is quadratic in nature. Our approach was to exploit the underlying group 
structure. Specifically, the polar decomposition of the dynamic variable al- 
lowed us to reduce the nonlinear initial value problem to a pair of linear-type 
evolution problems. Naturally, Theorem 3.1 implicitly makes a reference to 
Quantum Mechanics (via the operator W). In fact, the inter-connectedness 
of Quantum Mechanics and the Mesoscopic Mechanics will come to sharper 
focus in the next section. 
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4 General evolving- domain solution of the meso- 
scopic equation 

In the previous sections we have worked under the assumption that the do- 
main and image of the operator K defined in (1) are frozen in time. However, 
this assumption is neither necessary nor natural in the context of the meso- 
scopic equation (2). Indeed, it is natural to consider a more general setting 
when a priori both the domain and the image of operator K are allowed to 
evolve, i.e. 

K FG (t) : F(t) - G(t) (22) 

Here, it is understood that 

F(t) C Hi, and G(t) C H 2 

are finite-dimensional subspaces in two (possibly different) infinite-dimensional 
(separable) Hilbert spaces. In particular, the spaces F(t) and G(t) all inherit 
the Hermitian structure from the ambient Hilbert spaces. Furthermore, in 
this context, consider the Hamiltonian 

H(t) : .D — > Hi, 

which is well defined on a (fixed in time) dense linear subspace 

D C Hi. 

H(t) are also (formally) self-adjoint, i.e. 

(<p\H(t)ip) = (H(t)<p\ip) for all tp,ip £ D. (23) 

For a reason that will soon become clear we require a priori that 

F(t) C D (24) 

throughout the evolution. Finally, let us emphasize that the particular real- 
ization of the Hilbert spaces H 1 and H 2 as well as the Hamiltonian H(t) will 
remain implicit throughout our discussion as it is of no consequence to the 
conclusions we wish to draw. 

It is now clear how to interpret the mesoscopic equation (2) within this 
framework. Specifically, one needs to extend the operators K FG (t) through 
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zero to the orthogonal complement of F{t). Also, all operators K* FG and K FG 
need to be extended in an analogous way. Introduce the following shorthand 
notation 

K = K FG ®0 F ±, (25) 



and 

One checks directly that 
and, moreover, 



K~ l = Kp G ®0 G ±. (26) 



K' = K* FG ®0 G ±, (27) 



(K~ 1 )* = {K*)-\ (28) 



The shorthand notation seems intuitive and self-explanatory, and should not 
be confusing. We will refer to time-dependent families of operators of this 
type as the moving-domain operators. This terminology makes no reference 
to the 'moving image' as indeed, we will show that the image remains fixed 
for solutions of the mesoscopic equation, cf. Theorem 4.2. 

Definition 4.1 We say that a moving- domain operator K(t) as above is a 
local solution of (2) if for all ip G D ; equation 

{iUk + KH + B 2 (K*)- X )\$) = (29) 

holds for allt within a certain interval, say, t G [0, e). Of course, we write 

ihK = -KH - B 2 (K*)-\ 

The first goal is to show that the mesoscopic equation (2) has the unique- 
ness property even in this setting. In order to demonstrate this, the approach 
developed in the previous section will be exploited again. First, observe that, 
in view of (26), (27), and (28), equation (29) implies that for all e G H 2 and 
all i) G D 

f -ihK* + H*K* + B 2 K^) e\ib) = 0. 



Hence, the a priori assumptions (23) and (24) allow us to conclude that 

ihK* = H*K* + B 2 K~ 1 . (30) 
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Of course, the latter equation is understood in the ordinary sense. (If this 
may at first seem puzzling, let us point out that K* sends all vectors from H 2 
into F(t) C D. Therefore the equation can be 'evaluated' on all vectors from 
H 2 , and so it is expected to hold therein as, in fact, it does.) Next, observe 
that for an arbitrary e G H 2 , 

ihf t (KK*)e = 

= ihKK*e + ihKK*e 

= (-KH - B 2 (K*Y l ) K*e + K(H*K* + B 2 K~ 1 )e 
= —KHK*e - B 2 e + KH*K*e + B 2 e 
= 0. 

Here, the last equality is justified by the a priori assumptions (23) and (24). 
In particular, it follows that since G(t) is the image of KK*, it cannot evolve 
in time, i.e. 

G(t) = G(0). 

Therefore, it is possible to represent solutions in the polar decomposition 
with the self-adjoint and positive definite radial part R : G — > G, which is 
time independent. Now, suppose contrary to our expectation that (2) admits 
two a priori different moving-domain solutions on the interval t G [0, e), say, 

K (t) = RU (t), 

and 

K 1 (t) = RU 1 (t), 

while initially 

U o (0) = £/i(0). 

Here, 

Uo(t), C/i(i) : F(t) — >■ G, 

and the conventional extension to the whole space is understood implicitly. 
A direct calculation shows that 

ihU 0:1 = -Uo^H - H B U 0:1 , 

where H B = B 2 R~ 2 . Now, observe 

ih£ t (U U*) = -U HU* - H B U U* + U H*U* + U U*H B 
= -H B U Q Ut + U Q UtH B , 
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where, again, cancellation of two terms is justified by (23) and (24). At this 
stage, the extension of operator U{ to the whole of H 2 plays no role. In fact, 
we can view Y = UqU{ : G — > G as being the finite dimensional unitary 
operator satisfying 

ihj t (Y) = -H B Y + YH B . 

This is an equation of the type considered in Section 3, cf. equation (12). 
We already know it has the uniqueness property. Therefore, Y(t) = Id is the 
unique solution of this equation with the initial condition Y(0) = Id. Thus, 
U (t)Ui(t)* = Id, i.e. U (t) = Ui(t) for all t. In summary, we have 

Theorem 4.1 If the mesoscopic equation (2) in the broader moving -domain 
interpretation (cf. Definition 4-1) has a local solution K(t) in the interval, 
say, t G [0, e), then such a solution is uniquely defined by the initial condition 
K = K(0). 

The uniqueness property of (2) in such a broad Hilbert-space interpretation 
is a beautiful fact, indeed. Its proof relies on the inherent structure of the 
equation. 

Having established uniqueness of solutions we are empowered to find out 
the general form of solutions. Indeed, all we need to do is display a solution 
general enough to satisfy an arbitrary initial condition. Then, the uniqueness 
property will assure that no other solutions have been overlooked. This being 
the case, it would suffice to guess solutions, as long as they would be general 
enough. In what follows, it is shown how the general form of solutions can 
be deduced. 

In order to shed some light on the nature of moving-domain solutions, con- 
sider first a simpler case when F(t) = span {\ip{t))}, andG(t) = span {| </?(£))}, 
i.e. both spaces remain one-dimensional. Let operator K be represented in 
the form 

K(t)=a{t)\<p{t))m)\, (31) 

where a is a complex-valued function of time. Substituting, we find that 
equation (2) is translated into the following relation 

ih (%><V| + a\<p)(1>\ + a\<p)(i;\) = -a\ip){i>\H - ^>}<# 
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This latter equation is consistent if and only if there exist complex 
functions of time ci(t) and C2(i) such that 

ih\ip) = ci(*)|y>), 

ih{i>\ = y>\Mt)-H), 

and hence 

B 2 (t) 

iha = -(ci(t) + c 2 (t))a ^. 

Of course, the general solution of (32) is given by 

|^)>=exp(-ljT t Cl (Odt , )b(0)>. 
Next, introduce a new variable (^'(£)|, which is defined as follows 

<^(t)|=exp {^j\ 2 {t')dt^m)V 
The benefit of this is that 

Redefine also a by setting 

a = aexp (-Lj*( Cl (lf) + c 2 {t'))di S j . 

Observe that in particular K can now be re-written in the form 

K(t)=a'(t)\<p(0))m)\- 
Moreover, substituting (38) in (34) yields 

a'* 

Furthermore, setting a' = r exp (i$) leads to 

B 2 

Ihr — hr&) exp (i<£>) = exp 

15 



Now, since the exponential factor cancels, the real and the imaginary parts 
of the equation can be separated. In conclusion 



1 ft 

r = r = const, $ = — ^ / B 2 (t')dt' + %. 

flTf] J 





We summarize the result as 



Proposition 4.1 Consider operators K(t) : F(t) — > G(t), where F(t) and 
G(t) are one-dimensional spaces for all t. Equation (2) admits solutions in 
this form if and only if the following two conditions hold: 

1. The target space G(t) = G(0) = span{ \ip(0))} remains frozen in time. 

2. There is a vector tp' e Hi satisfying the one-particle Shrddinger equa- 
tion 

ih±W(t)\ = -W(t)\H, (41) 
which spans the domain spaces, i.e. 



F(t) = span{W(t))}. 



When both conditions 1 and 2 are satisfied, then the general solution of (41) 
admits representation in the form 

K(t) = roe** exp ^ jT B\t')dt^ |v?(0)> (<//(*)!, (42) 

where r and $ are arbitrary real numbers. 

I would like to emphasize that in particular the problem of existence of 
solutions of the nonlinear equation (2) has been reduced to the existence 
property of the linear Schrodinger equation (41). Naturally, the existence 
result and other properties of the latter equation are well known, e.g. cf. 
[8]. Moreover, in view of this result, even the notion of regularity of operator 
solutions of (2) acquires a clear meaning. 

Next, let us consider the general case of iV-dimensional domain and image 
spaces, which we will refer to as the N x iV-dimensional case. First, let 

F(t) =span{|^„(t)) :n= 1...N}, 
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and 

G(t) = span {\(p m (t)) :m — l... N} . 
Furthermore, let A be a complex matrix 

[0"mn(t)]m,n=l...N i Q"mn [^4]mn- 

Let the dynamic variable be represented in the form 

K(t)=Y,Ornn(t)\<P m (tW n (t)\- (43) 

(Summation is always carried out over repeated indices.) Observe that in 
particular 

Initially, some progress is achieved by exploiting analogy with the lxl- 
dimensional case. Indeed, observe that, in the N x iV-dimensional case, the 
mesoscopic Schrodinger equation (2) is translated into the following relation 

i%Y, (a mn \ip m )(il) n \ + a mn \ip m ){i/; n \ + a mn \(p m ) (ip n \^ = 

(45) 

Just as we have seen it before, also here a simple linear consistency check will 
help draw far-reaching conclusions. First, observe that for the equation to 
hold there must exist complex functions of time c' km (t) and c^(t) such that 

ih\(p m ) = Y. c 'km( t )Wk) for all m (46) 

and 

ih{4> n \ + (i) n \H = c?Mi\ for a11 n. (47) 
Secondly, introduce matrices 

C{t) = [c' km (t)km=l...N, C"(t) = [C(t)}n,l=l...N. 

Substituting (46) and (47) into (45), one obtains 

ihA + C'(t)A + AC"(t) = -B 2 (t) (A*)- 1 . (48) 

Note that operator K(t) as in (43) satisfies (2) if and only if the three condi- 
tions (46), (47), and (48) are satisfied by the |</?fc(t))'s, (■?/>„(£) |'s and the A(t). 
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In order to draw further conclusions, one ought to make the following ob- 
servations. First, one may require without loss of generality that both bases 
\ifi m (t)) and {ip n \ remain unitary during the evolution. Indeed, the operator 
K(t) can be described in arbitrary bases of F(t) and G(t). Now, suppose the 
two bases are unitary, say, at t — 0. Equations (46) and (47) imply that the 
bases will remain unitary for all time if and only if 

C'{ty = C\t), and C"{ty = C"(t), (49) 

i.e. if these matrices are Hermitian. With this understood, denote by V(t) 
and T"(t) the uniquely defined unitary matrices, which solve the two initial 
value problems: 

iht' = c'(t)r', r'(o) = id, 

and 

iht" = r"c"(t), r"(o) = id. 

In particular, 

\<Pm(t)) = J2iT'(t)] km \Mty} for all m. (50) 

Next, define a new unitary collection of vectors {ip' n (t)\ as follows 

mt)\=i:r'(t)} nl (Mt)\ for all n. (51) 

Naturally, the collection {ip' n (t)\ provides new unitary bases for spaces F(t). 
It has been selected in such a way as to simplify equation (47). Indeed, a 
straightforward calculation shows that 

m Jt m = f ° r aU H - (52) 

Furthermore, set 

A'(t) = V(t)A(t)T"(t). (53) 
and observe that (48) implies 

ihA' — -B 2 (t) (A'*) -1 . (54) 

We have already learned how to solve equations as this one in Section 3. 
Indeed, applying Theorem 3.1 (with H = 0) one obtains 

A'(t) = jA'(0)A'(0)* exp Q J* B 2 (t')(A' ^(O)A'(O)*)- 1 ^ . (55) 



18 



In fact, the initial conditions imply A'(0) = A(0). Finally, observe that in 
view of (53), (50), and (51) 

K(t) =E[A] mn \(p m (t)){Mt)\ 

= E[T'(ty 1 A'T"(ty 1 ] mn \ Vm (t)}(Mt)\ 

= E[A'] m „|^ m (0)>«(*)|. 

In the end, one ought to substitute (55) into the expression above. In sum- 
mary, we have the following result 

Theorem 4.2 Consider an operator K(t) : F(t) — > G(t), where F(t) and 
G(t) are evolving N -dimensional spaces. For K(t) to satisfy the mesoscopic 
Schrddinger equation (2) in the sense of Definition 4-1, it is necessary and 
sufficient that the following two conditions be satisfied: 

1. The target space 

G(t) = G(0) = span{\ip m (0)) : m = 1...N} 
remains frozen in time. 

2. There exists a collection {ip' n (t)} n=1 N which provides unitary bases for 
the domain spaces 

F(t) = span{\if>' n (t)) : n = 1 . . . N} . 

and, moreover, all the vectors satisfy the one-particle Shrodinger equa- 
tion, i.e. 

i% Jt m = ~ mH - (56) 

When both conditions 1 and 2 are satisfied, then solutions of the mesoscopic 
equation (2) admit representation in the form 

m=Y\ A \t)uwrnmmt)i (5?) 

where 

A\t) = ^'(0)^(0)* exp Q^ J B 2 (O(A / (0)A , (0)*)- 1 ^) . (58) 
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Naturally, Proposition 4.1 is a special case of Theorem 4.2. We introduced 
it beforehand not only because it is interesting in itself, but also because 
it provides a smooth introduction into the internal logic of the problem. 
As before, all strictly analytic issues, such as existence and regularity of 
solutions of the nonlinear problem (2) rest on the corresponding properties 
of the linear Schrodinger equation (56). Needless to say, vast literature is 
available in relation to the latter theme. 

I would also like to point out that Theorem 4.2 conforms with Theorem 
3.1. Both theorems expose a rather unobvious fact that the mesoscopic 
evolution is factored through the Schrodinger mechanics. Indeed, the W 
factor of Theorem 3.1 encodes the Schrodinger evolution. Of course, Theorem 
4.2 is not a mere corollary, and its proof required additional arguments, 
while also relying on the former theorem. Finally, let me emphasize that no 
conclusions have been drawn here as to the infinite-dimensional case, i.e. the 
case when (the nontrivial part of) the domain of K{t) cannot be encapsulated 
in a finite-dimensional space F(t). There are many other questions of interest 
not even attempted here, particularly those pertaining to the important case 
when the Schrodinger operator H itself depends on the dynamic variable K. 
The closing section explains the nature and significance of such a feedback. 

5 The broader context 

We will devote these concluding remarks to sketching the broader perspective 
in which the results of this article ought to be seen. First, I would like to 
point out that these results do not generalize to other types of operator equa- 
tions. Indeed, the unique type of nonlinearity in the mesoscopic Schrodinger 
equation plays a crucial role in the proofs. Specifically, it allows separation 
of the radial and the unitary part in the polar decomposition of the dynamic 
variable. Secondly, I would like to point out that the result is important in 
view of the physical interpretation of the MeM. Indeed, it shows that only 
the Schrodinger particles, i.e. electrons, can participate in the mesoscopic 
transport. This is not guarantied a priori, e.g. some other type of, say, a 
nonlinear wave could appear in place of the Schrodinger waves, which possi- 
bility is hereby a posteriori excluded. I will now provide a short synopsis of 
the physical interpretation of the solutions according to the theory that has 
been put forward in [5]. 
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Electromagnetic phenomena in vacuum are described by the classical 
Maxwell equations. These equations are modified by the so-called mate- 
rial constants or even by introduction of nonlinearities as modelers adapt 
them to describe propagation of the electromagnetic wave in various materi- 
als. Such an approach is usually sufficient when the model is meant to reflect 
what happens at the macroscopic scale. We know from experiment that in 
low temperatures some materials feedback to the electromagnetic field in a 
more profound way. Namely, at the nano-scale the spacial distribution of 
the magnetic field depends on the quantum picture of the electronic struc- 
ture of the material. This fact is of particular importance in the context of 
high-temperature superconductivity and the Quantum Hall Effects. A prob- 
lem arises, how to describe the inter-relation of the ambient magnetic field 
and the electronic structure. The particular form of this interrelation has 
far reaching consequences as regards the resulting galvanomagnetic charac- 
teristics of the material, e.g. cf. reference [2] in which some aspects of this 
problem are analyzed assuming random distribution of the magnetic field. 

The Mesoscopic Mechanics postulates that the distribution of the mag- 
netic field assumes a particular form depending on the quantum characteristic 
of the material. Specifically, let us focus attention on an idealized planar elec- 
tronic system characterized by the single-particle Hamiltonian H . Imagine 
this system being exposed to a perpendicular magnetic field with magnetic 
induction B. As a result of the interaction between the ambient magnetic 
field and the electronic structure, the magnetic flux will get distributed over 
the surface area nonuniformly. Here, depending on the properties of the 
system, the single-particle Hamiltonian H may or may not depend on the 
magnetic field. It is a basic precept of the MeM that even if the Hamilto- 
nian does not depend on the magnetic field, the flux distribution may still be 
nonuniform. Specifically, the MeM postulates that the distribution of flux is 
determined in a certain way by an operator K of the type considered in this 
article. Namely, let $ be the total magnetic flux through the surface, and 
let T stand for the coherent state 

T = E 

filled states 

which in this way accounts for the actual distribution of electron states. The 
simplest postulate of the MeM is that (with an appropriate normalization) 
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the distribution of the magnetic flux is approximately 

(x,t)^<S>\K(t)r\ 2 (x). 

Of course, the interpretation is probabilistic. Moreover, the evolution of the 
system is described via K by the mesoscopic Schrodinger equation (2). This 
equation of motion is determined by the total Hamiltonian S, cf. equation 
(3). The total Hamiltonian accounts for the single-electron portion of the en- 
ergy, as well as the inter-electron phase correlation energy. This latter energy 
is enclosed in the determinant (or entropy) term of S, which is switched on 
with an application of the magnetic field. I emphasize that we are looking at a 
new type of interaction of the magnetic field with the Fermi sea, independent 
and separate from the phenomenon of formation of Landau states. Electrons 
respond collectively since they are bound together by the energy of phase 
correlation. Naturally, the specific features of this phenomenon strongly de- 
pend on the energy-band structure of the material. Further information can 
be found in [5]. 

Finally, let me point out that the total Hamiltonian S is related to the 
following functional (whose arguments are functions) 

L A l$) = J |V^| 2 + 5 2 jln(H 2 ). (59) 

Naturally, the logarithmic integral mimics the entropy term. Some proper- 
ties of a particular realization of this latter functional, especially as regards 
magnetic- vortex type critical points, are described in an earlier article [7]. 
Last but not least, the MeM has a field-theoretic counterpart, the Nonlinear 
Maxwell Theory, cf. [6], which embraces (59) as one of its central objects, 
and provides models for many low-temperature phenomena. 
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